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Strange particle production at RHIC in a single-freeze-out modelB 


Wojciech Broniowski and Wojciech Florkowski 
The H. Niewodniczahski Institute of Nuclear Physics, PL-31342 Cracow, Poland 

Strange particle ratios and p± -spectra are calculated in a thermal model with single freeze-out, 
previously used successfully to describe non-strange particle production at RHIC. The model and 
the recently released data for 0, A, A and K* (892)° are in very satisfactory agreement, showing 
that the thermal approach can be used to describe the strangeness production at RHIC. 
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I. INTRODUCTION 


II. DEFINITION OF THE MODEL 


In this paper we calculate the ratios and the p_L-spectra 
for the strange particles produced at RHIC at yfs = 130 
GeV A. We use the thermal model with single freeze- 
out. The work is a direct follow-up of our previous study 
reported in Ref. [F]. The very recent data from the STAR 
collaboration on the production of <f> [§j, A , A [|| , and 
iC*(892)° [Q are confronted with the model and a very 
good agreement is found. 

Enhanced strangeness production in the high-energy 
nuclear collisions, relative to more elementary pp or e + e - 
collisions, was proposed many years ago j|U as a sig¬ 
nal of the quark-gluon plasma formation (for more re¬ 
cent arguments and discussion see, e.g ., [7]). In the 
meantime, the experimental evidence for the enhance¬ 
ment has indeed been found, however its significance for 
the plasma formation remains an open problem Q. The 
ratios of the particle abundances measured in Pb+Pb 
collisions at SPS (y/s =17 GeV A) and in Au+Au colli¬ 
sions at RHIC (y/s =130 GeV A) are very well explained 
in the framework of the thermal models which use the 
hadronic degrees of freedom only § 0 - In this approach 
one assumes that the net strangeness is zero, and the 
strange hadrons are in a complete thermal and chemi¬ 
cal equilibrium with other hadrons. A combination of 
the thermal model with a suitable hydrodynamic expan¬ 
sion (the single-freeze-out model of Ref. Q) led to a very 
good quantitative description of the RHIC transverse- 
momentum spectra of pions, kaons, and protons. In this 
paper we present the predictions of this model for the 
production of other particles, with a special emphasis on 
the strangeness production. Without any refitting of the 
two thermal and two expansion parameters of Ref. [|I|, 
we calculate the p_L-spectra of </>, A, A, iT*(892)°, S, E, 
and U. The spectra of </>, A, A and iF*(892)° are in 
a very good agreement with the very recently released 
data g-[|]. The spectra of other particles are also pre¬ 
sented and will be confronted with the incoming future 
data. We also compute the thermal spectrum of the J/if. 


The main assumption of our model is that the chemical 
freeze-out occurs simultaneously with the thermal freeze- 
out, i.e., the hadrons decouple completely when the ther¬ 
modynamic parameters reach the freeze-out conditions. 
In other words, the possible elastic rescattering processes 
after the chemical freeze-out are neglected. This assump¬ 
tion opposes the most popular scenario [Q, where the 
two freeze-outs are separated. The argument used in 
this context is the difference of scales, following from the 
fact that the elastic cross sections are, for most channels, 
much larger than the inelastic cross sections (note, how¬ 
ever, that this is not the case for pp interactions, where 
the inelastic channel dominates). The scenario with a 
single freeze-out is natural if the hadronization occurs 
in such conditions that neither elastic nor inelastic pro¬ 
cesses are effective. An example here is the picture of the 
supercooled plasma of Ref. Q . 

Recently, new hints have been presented in favor of 
smaller rescattering effects at RHIC: a successful re¬ 
construction of the iT*(892)° states has been achieved 
by the STAR Collaboration j|. As indicated by the 
authors, either the daughter particles from the decay 
iT*(892)° —> ttK do not rescatter or the expansion time 
between the chemical and thermal freeze-out is shorter 
than the iF*(892)° lifetime (r = 4 fm/c). In addition, the 
fact that the experimentally measured yields of iF*(892)° 
[|| are very well reproduced in the framework of the ther¬ 
mal models 00 suggests the picture with a very short 
expansion time between the two freeze-outs, as proposed 
in Ref. |]]. 

An important feature of our analysis is a complete 
treatment of the hadronic resonances, with all particles 
from the Particle Data Table [Q taken into account in 
the analysis of both the ratios and the spectra. This 
effectively leads to “cooling” of the hadronic spectra by 
30-40MeV [|ll| , which is a very important effect [^5| , cru¬ 
cial to obtain agreement with the data. 

The final ingredient of our model is the choice of 
the freeze-out hypersurface, which is, in the spirit of 
Ref. [p~G| [rif , defined by the condition 


r = Wf 2 - r/ L ~ r y ~ r z = const. (1) 
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The transverse size, p = y r% + r^, is limited by the con¬ 
dition p < p max . The £ and r z coordinates, appearing in 
the boost-invariant combination, are not limited, hence 
the model is boost-invariant. We have checked numeri¬ 
cally that this approximation works very well for calcu¬ 
lations in the central-rapidity region. Finally, we assume 
that the four-velocity of the hydrodynamic expansion at 
freeze-out is proportional to the coordinate (Hubble-like 
expansion), 



(i !jl H) 
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The question arises to what extent the assumptions 
(|l],p) are realistic. Typically, in a hydrodynamic ap¬ 
proach the freeze-out hyper surface contains, in the p-time 
plane, a time-like part, and a space-like part [|17| p3|. 
There is a conceptual problem here |2d| [2^j . A part of 
the particles emitted from the space-like surface goes 
backwards into the firecylinder and re-equilibrates. A 
common prescription is to exclude these contributions by 
hand. In our parametrization we neglect the space-like 
part altogether, thus avoiding the above problem. The 
time-like part of the hypersurface has, in many calcula¬ 
tions, the feature that the outer regions in the transverse 
direction freeze earlier than the inner regions. This is 
opposite to what follows from Eq. (1), or from the com¬ 
monly used versions of the blast-wave model, where the 
freeze-out occurs at a constant value of time in the p- 
time plane. Our present form for the hypersurface and 
for the velocity field corresponds to the so called scal¬ 
ing solution |l7],|2^,|2i| , which is obtained in the case of 
a small sound velocity in the medium. Naturally, valid¬ 
ity of the assumptions and their relevance for the results 
should be examined in a greater detail. In particular, the 
consequences of the specific choice of the freeze-out hy¬ 
persurface must be studied. The spherically symmetric 
case without the decays of resonances was investigated 
in Ref. |3(J. In Ref. [Q we have checked that two dif¬ 
ferent models lead to very close predictions for the spec¬ 
tra. Other parametrizations may be also verified with the 
help of the formulas given in the Appendix. One should 
bare in mind that at the moment there is no microscopic 
approach capable of reproducing all the features of the 
RHIC data (abundances, spectra, and HBT radii), such 
that the issue of how the hadronization and freeze-out 
occur is very much open. 

The model has four parameters possessing clear phys¬ 
ical interpretation. The first two parameters are the 
temperature, T, and the baryon chemical potential, p B . 
These are fixed by the analysis of the ratios of the particle 
abundances jn|]. The next two parameters, concerning 
the model of expansion, are the invariant time, r, and 
the transverse size, p max . The invariant time describes 
the lifetime of the system, and r 3 controls the overall 
normalization of the hadron multiplicities. On the other 


hand, the ratio p max /r influences the slopes of the p±- 
spectra. The expansion parameters have been fitted to 
the spectra IIIH in Ref. 0 by the least-square method. 
We note that p max is directly related, through Eq. (||), to 
the amount of the transverse flow, crucial for the shape 
of the transverse-momentum distributions. The maxi¬ 
mum value of the transverse fluid velocity, reached at 
the boundary, is /3™ ax = 0.66, while the average value is 
(P±) = 0.49. 

We use the previously-found values Q for the thermal 
parameters: 

T = 165 MeV, p B = 41 MeV. (3) 

The fit to the p± spectra of pions, kaons and protons for 
the most central collisions yielded [Q 

r = 7.66 fm, p max = 6.69 fm. (4) 

In the present paper we use the values of the parame¬ 
ters given by Eqs. (||) and (ft)), and calculate the spectra 
of other hadrons. As in Ref. [|j, the spectra are obtained 
from the Cooper-Frye [^8|,|4| formula 

/.&>•«), (5) 

where p^ is the four-momentum of the particle, is the 

volume element of the hypersurface defined by condition 
([]]), and fi is the phase-space distribution function for 
particle species i. It is composed from the initial and 
secondary particles, see the Appendix, proving that this 
is a correct approach for the case where ~ u^. The 
formulas valid for a general expansion are also given in 
the Appendix. 


III. RESULTS 

The thermodynamic parameters used in our model, 
Eq. ([|), yield the following ratios of the strange hadron 
abundances: 


n~ 

E~ £" 

(6) 


= 0.18, — = 0.55, —= 0.20, 

Zj YY 


= 0.47 (0.49 ± 0.03) exp t.) 

(7) 
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+ 
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= 0.85, — = 0.76, — = 1.02. 
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The experimental value for the A /p ratio follows from the 
data published in Refs. p|,|33|. The feeding from weak de¬ 
cays has been included in the above ratios. For example, 
the multiplicity of S contains a contribution from the Q 
decay, and the multiplicity of A contains contributions 
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from the decays of £, X, and (indirectly) fl. When the 
feeding from weak decays is excluded, the following ratios 
are noticeably altered: 

y- ^ 

— = 0.36, - = 0.26. (9) 

A p w 

The abundance of the £ hyperon is sensitive to the inclu¬ 
sion of higher £ resonances, in particular £(1690), whose 
properties (spin) are not known. This leads to a theoret¬ 
ical error of at least 10%. This is a general feature of the 
thermal fits. The heavier particles receive contributions 
from their excited states, however the experimental par¬ 
ticle spectrum is less known as the mass increases. The 
higher mass states are suppressed by the thermodynamic 
factor, on the other hand, they are more numerous, ac¬ 
cording to the Hagedorn hypothesis. For the recent dis¬ 
cussion of this issue see Refs. [ ]35| , |36| . 

For the J/'ip we obtain 

= 0.12 x 1CT 6 . (10) 

7r 

The abundance of J/ijj is an order of magnitude lower 
than the values obtained in Refs. |37],|3^], where a wide 
class of reactions is studied and the ratio (J/^)/h~ ~ 
10 -6 is found. Note, however, that due to the large mass 
of J /^5 the ratio is very sensitive to the chosen freeze- 
out temperature. The ratios dHH) are predictions of 
our model, to be verified by future data. More elaborate 
models, taking into account non-thermal J/'ijj produc¬ 
tion, are described in Refs. f39[ , [40|j . 

More detailed information about the hadron produc¬ 
tion is, of course, contained in the transverse-momentum 
spectra. Several calculations have been performed re¬ 
cently in the hydrodynamic approaches mm- In the 
framework of the thermal models such a calculation is 
much more involved than the analysis of the abundances, 
since it requires the implementation 00 of the cascade 
processes according to the formulas of the Appendix. The 
results are shown in Fig. [l|. The upper part displays the 
spectra of pions, kaons, antiprotons (already presented 
in 00), the </> mesons, and the IT*(892)° mesons. With 
the parameters fixed in Ref. Q with help of the spectra 
of pions, kaons, and antiprotons, we now obtain the p±- 
spectrum of the </> mesons, which agrees very well with 
the recently-reported measurement [|| . The model curve 
crosses five out of the nine data points; one should also 


Tn Ref. JI| the preliminary data for p were presented 
that exclude the feeding from the weak decays through the 
use of the HIJING model. The correction for weak decays 
results in about 20% reduction of the normalization of the 
spectrum. In the present paper the official data are presented 
[p3|, which include the full feeding from weak decays. 


bear in mind that the systematic experimental errors are 
expected at the level of about 20% 0. The experimental 
ratios for the most-central events [2] agree well with the 
output of our model: 

£=0.019 (0.021 ±0.001)expt., 

T 73 = 0.15 (0.10 - 0.16)expt, (11) 



FIG. 1. (a) The p^-spectra at midrapidity of pions, kaons, 

antiprotons, the 0-mesons, and the K*( 892)° mesons, and (b) 
of A, A, X - , £ _ , Q - , and J/'i/j. The open and filled sym¬ 
bols describe the STAR and PHENIX highest-centrality data, 
respectively, for Au + Au at s 1//2 = 130 GeV A. All theoret¬ 
ical curves and the data are absolutely normalized. The two 
expansion parameters of the model were fitted to the spec¬ 
tra of pions, kaons, and antiprotons in Ref. [1]. The model 
curves include full feeding from the weak decays. The data 
come from Refs. [2,3,31-33]. The STAR data for Tr~,K~,(j) 
and AT*’s are preliminary. The data include full feeding from 
weak decays. 

We recall that the </> meson deserves a particular at¬ 
tention in relativistic heavy-ion collisions. It serves as a 
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very good “thermometer” of the system, since its interac¬ 
tion with the hadronic environment is negligible. Also, it 
does not obtain any contribution from the resonance de¬ 
cays, thus its spectrum reflects directly the distribution 
at freeze-out and the flow. 

The upper part of Fig. |l] shows also the averaged spec¬ 
trum of K*' s. In this case the model calculation is com¬ 
pared to the very recently released data ||. Once again 
we observe a very good agreement between the model 
curve and the experimental points. The measured ratios 
involving the yield of K* (892)° [Q also agree well with 
the output of the thermal model: 


0.90 (0.92 dz 0.14) expt ., 

0.046 (0.042 ± 0.004 db 0.01) exp t., 

0.33 (0.26 ± 0.03 ± 0.07) exp t., 

0.42 (0.49 db 0.05 ±0.12) expt .. (12) 

As already mentioned in the Introduction, the success¬ 
ful description of both the yield and the spectrum of 
if* (892)° mesons supports the concept of the thermal 
description of hadron production at RHIC, and brings 
evidence for small expansion time between chemical and 
thermal freeze-outs. If the if* (892)° mesons decayed be¬ 
tween chemical and thermal freeze-out, the emitted pions 
and kaons would rescatter and the if* (892)° states could 
not be reconstructed. In addition, if only a fraction of 
the if* (892)° yield was reconstructed, it would not agree 
with the outcome of the thermal analysis which provides 
the particle yields at the chemical freeze-out. Thus, the 
expansion time between chemical and thermal freeze-out 
should be smaller than the if* (892)° lifetime, r = 4 fm/c 

0 . 

The lower part of Fig. |lj shows the spectra of strange 
baryons, and the J/'ip meson. The data for A and A were 
taken from Ref. [||. The p ±-dependence of the spectrum 
of A and A is quite satisfactorily reproduced, within 20%, 
by the model calculation. The other curves in the figure 
are our predictions. 

In Fig. |]we give the inverse slope parameters for var¬ 
ious particles, defined as 


o 


A; = — 


■ In 


dNi 


dm_\_ \ 2iTm_\_dm_\_dy 


y= 



if* 0 

if* 0 

h~ 

if* 0 

~K~ 

if* 0 


practiced in the interpretation of the data. Secondly, 
the inverse slope increases with the growing mass of the 
hadron, however this increase is not linear and it is not 
possible to parameterize it with a simple formula. Re¬ 
call that the spectra are fed by the decays of resonances, 
which is an important effect, significantly reducing the 
inverse slope at lower values of m± [|lT| . The m±_ spectra 
of pions and kaons are convex [j], whereas for the remain¬ 
ing particles they are concave, the more the heavier the 
particle is. All curves in Fig. asymptote to the value 

Aoo ~ , = -= 363MeV. (14) 

V 1 + Prna x/t 2 - Pmax/T 

This formula, specific for our expansion model, is simple 
to obtain from Eq. (^j) with fi given by the exponential 
form, which is true at large m_ j_, where the feeding from 
the resonance decays is absent. As can be seen from 
Fig. |, the asymptotics is reached slowly for most of the 
particles. 



FIG. 2. The inverse-slope parameters for various particles, 
calculated as a function of the transverse mass, m±. The 
constant mo is the mass of the given particle. 


IV. EXCLUDED-VOLUME EFFECTS 


This definition reduces to a constant for the exponential 
spectrum, exp(— m±/X). As can be seen from the figure, 
the inverse slopes depend strongly on the value of m± 
in the region of interest. The same conclusion has been 
reached in Ref. |20|. Thus, at least from the viewpoint 
of our model, it is impossible to describe the spectra in 
terms of a single slope parameter, as is most commonly 


Our fitted values for the geometric parameters r and 
Pm ax are low, of the order of the size of the colliding nu- 


i We use the convention that a function with positive second 
derivative is convex. 
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clei. This leads to two problems: 1) the values of the 
HBT radii |I^] are too small compared to experiment, 
and 2) there is little time left for the system to develop 
large transverse flow. The problems can be helped by the 
inclusion of the excluded-volume (van der Waals) correc¬ 
tions. Such effects were realized to be important in the 
earlier studies of the particle multiplicities in relativis¬ 
tic heavy-ion collisions p3^48[, where they led to a sig¬ 
nificant dilution of system. In the case of the classical 
(Boltzmann) statistics, which is a very good approxima¬ 
tion for our syste m p9| , the excluded volume corrections 
bring in a factor [47|] 


e -Pvi/T 

1 + Vje~ Pv i/ T rij ’ 


(15) 


into the phase-space integrals, where P is the pressure, 
Vi = 4|7rrf is the excluded volume for the particle of 
species i, and rii is the density of particles of species 
i. The pressure is calculated self-consistently from the 
equation 


P = E ^°( T ’ W - Pv i/ T ) = E ^°( T ’ Mi)e- P " i/T , 

(16) 

where Pf is the partial pressure of the ideal gas of 
hadrons of species i. For the simplest case where the ex¬ 
cluded volumes for all particles are equal, ri = r, the cor¬ 
rection manifests itself as a common scale factor, which 
we denote by S~ 3 . The Frye-Cooper formula can then 
be written in the form [Q 


dNj 

d 2 p_\_dy 


/ + °° rPmax /T 

da\\ / sinha^d (sinhajJ 

-oo J 0 


nZ 7T 

x / d£p ■ uS~ 3 fi(p ■ u) , (17) 

Jo 

where p • u = mxcosha||Cosh<a_L — p_ l cos £ sinh<a_L. The 
emergence of the factor S~ 3 in Eq. (17) may be compen¬ 
sated by rescaling p and r by the factor S. That way the 
system becomes more dilute and larger in such a way, 
that the particle multiplicities and the spectra are left 
intact. 

For our values of the thermodynamic parameters, with 
pi) = 80MeV/fm 3 , we find S = 1.3 with 
r = 0.6fm, and S = 1.6 with r = 0.8fm. Such values of 
the excluded volumes have been typically used in other 
calculations. That way, the increase of the size param¬ 
eters at freeze-out of the order of 30-60% is generated. 
Therefore, the inclusion of the excluded-volume correc¬ 
tions opens the possibility that the problems 1) and 2) 
may be alleviated. More detailed analysis will be pre¬ 
sented elsewhere. 


V. CONCLUSIONS 

The presented model results for the strange hadron 
production at RHIC support the idea that also these par¬ 
ticles are produced thermally. No extra parameters (e.g. 
strangeness suppression factors) are necessary. This is an 
important information concerning the particle produc¬ 
tion mechanism at RHIC. Moreover, the hypothesis of a 
single freeze-out of Ref. U has been additionally verified 
with the available spectra of 0, A, A, and if* 5 s. In short, 
the thermal model supplied with a suitable parameteriza¬ 
tion of the freeze-out hypersurface and velocity (i.e. the 
inclusion of the longitudinal and transverse flow) allows 
for an efficient and uniform description of the RHIC spec¬ 
tra. A further verification of the model will be provided 
by the calculation of the HBT correlation radii, as well 
as the analysis of non-central events, which are currently 
being examined. 

We are grateful to Marek Gazdzicki for helpful discus¬ 
sions. 


APPENDIX A: 


Let us consider a sequence of the resonance decays. 
The initial resonance decouples on the freeze-out hy¬ 
persurface at the space-time point xjy, and decays af¬ 
ter an average time inversely proportional to the width, 
tat ~ 1/Tat- We follow one of the decay products, formed 
at the point xn-i- It decays again after a time tat-i, and 
so on. At the end of the cascade a particle with the label 1 
is formed, which is directly observed in the experiment. 
The Lorentz-invariant phase-space density of the mea¬ 
sured particles is (we generalize below the formula from 
Ref. @ where a single resonance is taken into account) 

ni (xi,pi) = 


/ 



dr 2 r 2 e ~ r2T2 


. P2T2 

x 2 H- Xi 

m 2 


J B (Pn,Pn- i) J drN^Ne Tntn 


x J (x N ) p^ N £ (4) ^xat 
x In [pn • u (xat)] , 


Pnt n 

m N 


- XN-l 


(Al) 


where B(q, k ) is the probability for a resonance with 
momentum k to produce a particle with momentum g, 
namely 


B(q,k) 



m R 


(A2) 
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Here b is the branching ratio for the particular decay 
channel, and p* (E*) is the momentum (energy) of the 
emitted particle in the resonance’s rest frame. Integra¬ 
tion over all space-time positions gives the momentum 
distribution 

„ dN x [ 4 

Epi d?pl = J d XlUl 

= f -jJ—B (p 2 ,Pi) ■■■ f -77—-B (p N ,p N ~i) 

J ^P2 J ^PN 

x / ( x N ) p p N f N [p N • u (xat)] , 

(A3) 

which should be used in a general case. Note that the 
dependence on the widths Tk has disappeared, reflect¬ 
ing the fact that it is not relevant when or where the 
resonances decay. 

A simplification follows if the element of the freeze-out 
hypersurface is proportional to the four-velocity (this is 
exactly the case considered in our model defined by Eqs. 
( 0 ) and (§)), 


dT, fl (xN) = <1Ti{xn) u^xn)- (A4) 


Then 


E„ 


(IN i 
' Pl Wi 
d 3 p N 


= JdN(x N ) j ^*B(p2,pi) 


B (pn,Pn- i) Pn • u (xn) /tv [pn • u (xn)} 


*/ 

= y dS ( x N ) Pi • u (x N ) fi[pi-u (x N )} , (A5) 

where we have introduced the transformation 


Pk-l ■ u (Xjv) fk -1 \pk -1 • U (zjv)] 
d 3 Pk 


-I 


Eri 


-B (pk,Pk-l)Pk • u ( XN ) fk [Pk • u (xn)} , 


(A6) 


which can be used (step by step along the cascade) to cal¬ 
culate the distribution of the measured particles. In the 
fluid local-rest-frame, most convenient in the numerical 
calculation, we have #(%) = (1,0, 0,0), and the trans¬ 
formation ( |A6| ) reduces to the form discussed in |ll|. A 
technical simplification relies in the fact that in Eq. ( [Aq) 
the space-time integration over the hypersurface is per¬ 
formed at the end, consequently the momentum integra¬ 
tion in Eq. (A6) preserves the full spherical symmetry. 
That way, the momentum integrals are one-dimensional 
and the numerical procedure is very fast. On the other 
hand, in the general case of Eq. (A3) the integration 
over the hyper-surface has to be done first, reducing the 
symmetry of the following momenta integrals. The sym¬ 
metry is cylindrical, thus the momentum integrals are 


two-dimensional. Moreover the integrands have an inte¬ 
grate singularity, which make the numerical procedure 
more involved 00- 


* Supported in part by the Polish State Committee for 
Scientific Research, grant 2 P03B 09419. 
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